Efficient control of transient wave forms to 
prevent spreading depolarizations 
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Abstract 

In various neurological disorders spatio-temporal excitation patterns constitute ex- 
amples of excitable behavior emerging from pathological pathways. During migraine, 
seizure, and stroke an initially localized pathological state can temporarily spread in- 
dicating a transition from non-excitable to excitable behavior. We investigate these 
transient wave forms in the generic FitzHugh-Nagumo (FHN) system of excitable 
media. Our goal is to define an efficient control minimizing the volume of invaded 
tissue. The general point of such a therapeutic optimization is how to apply control 
theory in the framework of structures in differential geometry by regarding param- 
eter plane M of the FHN system as a differential manifold endowed with a metric. 
We suggest to equip M with a metric given by pharmacokinetic-pharmacodynamic 
models of drug receptor interaction. 
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1 Introduction 



Excit ability is a central concept in neurophysiology both in health and dis- 
ease ( iHilld . 119971 ). Under some pathological conditions altered excitability can 
spread through the tissue. Such tran s ient wave forms (TWF) occur for exam - 

ple in mi graine with aura (Lauritzen . 1994 : Aurora et al . 1998 : Welch . 2005[). 

epilepsy (Somien fc Aitken . 1984 Giorgi et a/. . 200611. and stroke (Nedergaard fc Astrup , 
19861 : iHossmannl . Il994l : iFisher fc Takanol . Il995l : iHakinil . Il998l : iBackl . Il998h 71n 
migraine with aura seizure-like activity spreads slowly through parts of the 
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cortex. This is observed by sym ptomatic (ILashleyl . Il94ll ) and electrophysio- 
logical (jHadjikhani et all 1200 ll ) events. Likewise, epileptic seizures can have 
a localized onset and then grow in intensity and start to spread. This usu- 
ally leads to subsequent generalized motor involvement commonly referred to 
as partial seizures with secondary generalization. In some cases, however, the 
epileptiform activity may induce changes to subcortical structures produc- 
i ng clinical signs of ge neral motor involvement that merely mimic a spread 



( ISchindler et all 120071 ). Last but not least, during stroke a cortical region 



that surrounds the infarct core and that initially suffers functional injury 
can gradually and pro gressively fail and suffer irreversible structural injury 
in untreated patients (IBaird et all 119971 ). In this condition the electrophysi- 
ological changes are persistent, but there exist intermediate repetitive TWF 
and it was suggested that 'therapy might [...] target the intermediate forms 
of spreading depolarizations so as to protect the penumbra [tis sue surround' 



ing t he infarct core] against recruitment into the infarct core' (jPreier et al. 
20061 ). These three paradigmatic clinical manifestations of spreading patholog- 



ical states motivate efforts to understand how the spread of such states arises 
and how it can be controlled. 



The phenotypical manifestation of such s preadi i ig pa thological states is the 
phenomenon of spreading depression (SD) (ILead . Il944l ). SD is the basis of mi- 
graine with aura. Essentially identical electrophysiological features are associ- 



ated with infarct expansion, called peri-infarct dep olarization (PID) (IMayevsky et al. 



19961 : IStriimpel et all l2002l : iFabricius et al[\200w. Th e cascade of events that 
produce SD is related to seizures (iKager et a/TEoOOh . Despite certain differ- 
ences, seizure events that begin with epileptiform discharges can either termi- 
nate in SD, facilitate the synchronization, or spread by a similar mechanism 



over a large area with a velocity resembling that of SD (iBures et all 11974 



Somjen fc Aitkeru . Il984j : iGorji &: Speckmannl . 120041 ). SD is seizure- like activ- 
ity evolving as a slowly spreading non-or-all type process. It is characterized 
by the feedback of ion currents that ch ange ion concentra t ions, which, in 



turn, influence the membrane potential (IKager et all l2000l : ISomjenl . 120011 ). 
Shortly after its onset all neuronal activity is depressed, hence its name. 
The name is misleading, because SD can still be observed even when neu- 
ronal activity is depressed b y blocking the fast transient sodium current l^a 
(ITobiasz fc Nicholsonl . Il982l ). SD emerges from an excitable pathway in neu- 
ronal tissue independently of the normal neuronal activity. It was therefore 
suggested to cat egorize SD a r id sim ilar phenomena under the term spreading 
depolarizations (IDreier et a/.l . l2006l ). 



There is ample evidence that SD belongs to the self-or ganization processes due 



to the coupling of biochemical re action with diffusion (IMartins-Ferreira et al. 



gested (Tuckwell fc Miural. 



Kager et all l200d : ISomjen 



1978 



2000l : iDahlem fc Chroni eld . 120041) . Mathematical models of SD have been sug 



Reggia fc Montgomery! . Il994j : IShapird . 12001 



200ll ). though there is not yet consent on the 
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mechanism. We will use the spatially extended FitzHugh-Nagumo (FHN) sys- 
tem as a generic model of neuronal excitation patterns based on reaction- 
diffusion. As a neuronal model, it describes generic pattern formation prop- 
erties not limited to nerve impulse propagation along an axon, although it 
was originally deriv e d fro r a the Ho dgkin- Huxley model of action potentials 
flHodgkin fc Huxlevl . Il952l : IPitzHughl . [l96ll : iNagumo et oZI . Il962h . T he FHN 



syste m describes also the spatial features of SD wave in neuronal tissue (jPahlem &: Miiller 
20041 ). Furthermore, the transition from non-excitable to excitable m edia sup- 
porting traveling wave s was well investigated in the FHN model (jWinfred . 
199ll : iHakim &: Karmal . Il999l ). It was suggested that th e spati o-temporal pat- 
terns in SD occur at this transition (jPahlem fc Miilleii . |2004| ). The regime in 
which this transition takes place is also well investi gjated in chemical model 
syste ms in experiment and theory, for a review see (iMikhailov fc Showalten . 
2006h . 



The route to spreading depolarizations in a generic model is provided by two 
independent pathways: one lowering the threshold of evoked pathological ac- 
tivity, the other changing the time scale of biochemical reaction rates, and 
hence their time scale ratio. The two pathways might offer new opportunities 
in developing optimal therapy. Consider the case that a pathological condi- 
tion is caused by a shift along one path whereas therapeutic strategies are 
available for both pathways. Suppose both strategies can be combined while 
each has individual response rates and side effects. What is an optimal time 
efficient combined therapy stopping the spread while minimizing side effects? 
Since therapies can be combined, there is a two dimensional manifold in which 
therapy takes place. The strategy we suggest is to equip this manifold with a 
metric that allows us to find an effective combined therapy with minimal side 
effects. Effectiveness is defined by finding a path to a sufficiently low excitabil- 
ity where the tissue is not susceptible to spreading events, whereas efficiency 
refers to side effects and time. 



2 Parameter space of the FHN system 



We assume that a standard activator-inhibitor scheme leads to the observed 
propagation phenomena during SD. The activator and inhibitor variables, u 
and V, are coupled by their kinetic reaction rates f{u,v) and g{u,v), respec- 
tively, and can diffuse in the medium. The equations are 
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Diffusion is represented by diffusion coefficients Du and D^. By re-scaling 
space, the ratio 5 = of the diffusion coefficients can be introduced replac- 
ing the parameters and D^. The parameter e is the time scale ratio of 
inhibitor and activator variables. The reaction rates f{u,v) and g{u,v) may 
possibly be derive d from a more complex model of SD, e.g., the one from 



Kager et all ( 120001 ). by lumping together all activator variables, such as in- 
ward currents and extracellular potassium concentration [-^''''jo into a single 
activator variable and their combined kinetics into a reaction rate /(■,■). Like- 
wise, a single inhibitor variable could be related to recovery processes, such 
as effective regulation of [-ft'^jo by the neuron's Na-K ion pump and the glia- 
endothelial system. This will be an important task of future investigations. 
We aim to describe universal features of reaction-diffusion coupling that lead 
to the onset of spreading pathological states and do not specify the variables 
u and V which underly these characteristics other than that they play the 
roles of activator and inhibitor, respectively. Their kinetic functions f{u,v) 
and g{u, v) are given by the FHN system 

f{u,v)=u-^u^ -V (2) 
g{u, v) = u + P — •yv 

where f3 and 7 determine the excitation threshold (Fig. 1). 

We shall start by viewing the parameter plane of the FHN model as a manifold 
M and consider its geometric structure. M has four dimensions (e, /3,7,(5). 
In a certain regime the parameters f3 and 7 determine the threshold for a 
non-or-all excitation process. Like e, their variation can cause a bifurcation: 
the emergence of sustained travelling waves. Instead of the four dimensions 
usually a two dimensional subse t is investigated to describe t his bifurcation, for 



exarn ple the section at 7 = 0.5 (jWinfred . Il99ll ). or at 7 = (IHakim fc Karma 



19991 ). both with 6 = 0. Consider the subset at 7 = 0. A particular FHN 
system is specified by a point q of this subset. It can be parameterized by the 
coordinate functions e(g), i.e., the time scale separation of u and v, and the 
threshold parameters /3(g). As an alternative coordinate function for (3{q) 

A{p) = {(3-1(3') (3) 



can be chosen (IHakim fc Karmal . Il999l ). While /3 is a measure of the threshold. 



A is related to a measure of excitability, because it is equal to the inhibitor 
concentration in the steady state (Fig. 1). This is rather a convention than 
a definition of excitability. To be more general, we shall only assume that 
excitability E is a C°°-function :[/—>■ M in a subset U of M. Further prop- 
erties of this function will be defined later. Firstly we want to note that there 
obviously exists a whole set of coordinate systems A for M with coordinate 
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transformations like Eq. (3) being C°° diffeomorphisms. Naturally, M can be 
identified as a C°° differential manifold. For a given subset U of dimension 
n we can choose the coordinate system ^ = f^"*^, . . . , = [C] ^ -A., with 
i G {1, 2, . . . , n} that seems to suit best the purpose of study. For example, in 
the section with n = 2 and = e the parameter can be chosen either as (3 or 
as A. We will consider the hypersection at 5 = 0, and 7 = and use various 
coordinate systems in this submanifold. 



3 Transient Wave Forms (TWF) 



Excitability is an emergent property of active media. It arises when a critical 
parameter value is crossed above which the medium is sus c eptibl e for sus- 



tained propagating excitation patterns (Mikhailov fc Zvkovl . Il99ll : IWinfree 



19911 : iHakim fc Karmal . Il999l : iMihahuk et all 120021) . In a ID medium this bor- 
der in M is the propagation boundary dP. It is obtained by finding solutions 
in a co-moving frame in the form u{x,t) = U{y), v{x,t) = V{y), z = x + ct, 
where c is an additional parameter of the unknown wave propagation speed. 
Travelling wave solutions in the form of a single pulse are thus equivalent to 
the existen ce of a hornoclini c orbit satisfying a system of ordinary differential 
equations (IKuznetsovl . 1 19951 . p. 224 ff). The transition marks a bifurcation of 
codimension one. Therefore dP is a hypersurface in M separating the regime 
supporting travelling waves from the non-excitable one. In a two-dimensional 
parameter space, spanned by e, and P or A, dP is a curve (Fig. 2). Below dP 
any confined perturbation of arbitrary profile decays eventually. Above dP 
some wave profiles are stable and travel with constant velocity. 

Consider the region close to dP, which is referred to as sub-excitable. There 
transient wave forms (TWF) exist though sustained waves may not. Note first 
that in the absence of a metric in M closeness to dP can only be defined with 
reference to other borders by which M is further subdivided. The adjacent 
bifurcation curve above dP is dR. It is only defined in systems with more than 
one spatial dimension. In these systems dR is the border above which open 
wave fronts will not disappear because each open end curls in to form a spiral. 
Spiral waves rotate and re-enter multiple times the invaded tissue. Below dR, 
but still above dP, TWF exist depending on the initial size o f the excitation . 



In 2D such waves occur in the shape of particle-like w aves (jMihaliuk et al. 



20021 ). which can be controlled by excitability gradients (jSakurai et a/.l . l2002l ). 
In this study, we investigate the region below dP where unstable TWF exist 
in the FHN systems with one spatial dimension. 

To further subdivide M in the regime below dP, we suggest to take the dis- 
tance a TWF spreads (Fig. 3b) to define isolines. This distance defines the 
volume of tissue at risk (TAR) referring to the risk of transient neurologi- 
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cal symptoms or even of permanent damage (PID case) when cortical tissue 
is invaded following a local stimulation. At dP stable traveling waves exist. 
They can invade the whole tissue. Therefore the TAR value is infinite at dP. 
Below dP the value of TAR is finite and rapidly decreases to zero when one 
moves the FHN system away from dP in parameter space (Fig. 3b-c). We 
found the existence of a new boundary dS as the TAR value approaches zero. 
At this boundary any stimulation profile collapses into the steady state with- 
out broadening, that is, without affecting surrounding tissue. dS is therefore 
contrary to dP at which the TAR value is infinite. The region between ^-S" 
and dP is the region where TWF exist in a ID system. This region defines a 
subset U of M. The general point of a therapeutic control strategy is to leave 
U by crossing dS. 



4 Efficient control of transient wave forms 



So far excitabihty E has not been defined. We assume £■ to be a C°°-function in 
U. It is reasonable to assume that E is constant on the two borders, because 
these are bifurcation lines at which TAR is constant. Furthermore, since a 
change in TAR is suggestive of a change in excitability, we propose both to be 
linearly related. 

Our goal is to control a path A in [/ along which excitability, and with it the 
volume of TAR, is efficiently diminished. Let the path A start mU ai q where 
the neuronal tissue temporarily supports the spread of excitation. Let the path 
end on dS. Firstly, we assume sufficient knowledge of control characteristics, 
e.g., by pharmacokinetic and pharmacodynamic means, as described in the 
next section. We will show in the next section that this provides a metric on 
U . In this section we will just assume U to be a differential manifold with 
a metric. Furthermore, we assume (i) that our control method allows us to 
choose any path A : 7 — > [/ parameterized by some interval / C M, (ii), as 
already stated, that £■ is a C°°-function in a subset U, which includes A, and 
(ni) that E{q) > E{dS) holds. 

Which path should we take, if we want to reduce excitability by going from q 
to dS by deliberate control? The efficiency of successful control dragging the 
system into the target state p on dS should be given by some optimization 
criterion. When a metric is given two paths are privileged: the shortest path A^ 
between q and dS^ and also the one that minimizes E by gradient descent. The 
latter path A^ implies a metric because covectors like dE/d^^ (representing a 
gradient) and contravariant vectors (represented by tangent vectors to a path 
A) are unrelated objects of different kinds. Only a metric tensor g^j defines 
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the gradient as a tangent vector (using summation convention) 



(grad^;)^ = g- 




(4) 



Therefore, only when a metric is given on [/, we can apply some optimization 
criterion for the efficient therapeutic path. 



5 Metric tensor on U 

As efficient control, we proposed in the last section a method based on a 
metric structure in U . Therefore, a metric structure is needed to optimize 
control. Hence it is natural in this context to endow U with a metric that is 
derived by some sort of cost function of the control method. We introduce 
a standard pharmacokinetic and pharmacodynamic scheme to illustrate this 
concept. Let Q be the concentrations of drugs which regulate diverse functions 
in populations of neurons. For the sake of simplicity, we neglect the details 
of pharmacokinetics as the discipline that describes dosage regimes and the 
time-course of Q in the body by absorption, distribution, metabolism, and 
excretation. We assume the drugs can be constantly administered and their 
rate of administrations equals their rate of metabolism and excretion. Thus 
Q is immediately in its steady state vahic. Furthermore, wc assume that Q 
follows linear pharmacokinetics. In this situation the steady state of changes 
proportionally according to dose. 

The relation between drug dose and response is usually modeled as a hy- 
perbolic function assuming a simple drug receptor interaction. Suppose the 
response to Q is 



where the define a new coordinate system ^ = [^^, . . . , in U. In this equa- 
tion, r* denotes transducer functions that represent the response of the FHN 
system to the drug Q. For the sake of simplicity we use as transducer functions 
r* the identity. ECIq are the effective doses 50, i.e., doses at which 50% of the 
maximal responses ^l^^^ are achieved (Fig. 4a). ^l^^^ is the asymptotic value 
of for large concentrations. 

In analogy with Eq. (3), where we have introduced the new coordinate A, we 
have thus introduced the coordinate system ^ G ^ as new control parameters 
of the FHN system. A new coordinate system with a known metric is needed 




(5) 
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as a reference system from which the metric in U is obtained by coordinate 
transformation. Let us denote the components of the metric tensor of U in 
the coordinate system C, as Qij. Above we have introduced the new coordinate 
system ^ G ^ of the response variables. The components of the metric tensor 
in this coordinate system (^) are then 



Likewise, we can calculate the unknown components gij from gij. The coordi- 
nate system is a privileged reference system for U by selecting However, 
this does not imply that it is Cartesian [gij = 6ij, 6ij being the Kronecker 
delta tensor representing the Euclidean metric in Cartesian coordinates), as is 
readily seen since the choice of concentration units is arbitrary. For example, 
commonly the logarithm of the concentration is plotted on the abscissa, 
as shown in Fig. 4a. Hence the question is whether a structure of the phar- 
macodynamic scheme can be natural choice treated as a metric structure in 
differential geometry defining g^. 

Suppose that C o^' oiis of their metabolites have toxic side effects. Their dose 
response curve follows the same relation as Eq. (5), although shifted to the 
right on the dose axis by a higher toxic dose 50 (TC50) 

• / C \ 

rpt It I max^ I 



f denotes transducer functions that represent the response of the toxic system. 
Again for the sake of simplicity we use as transducer functions f the identity 
function. In analogy to Eq. (6) we obtain 

_ _ dC dT^ (9T' dC^ 
9»P - d^~dC WW' 



assuming T = [T^, . . . ,T"] builds a Cartesian coordinate system of the costs 
with the Euclidean metric 6ki. The association of the cost coordinate system 
(T) with the Cartesian one is a naturally choice, since the efficiency of control 
is-by definition-measured in these costs. 

We end this section with one example that to some extent corresponds to a 
subtype of migraine with aura. For this subtype a novel pathogenic genetic 
mutation was found. Therefore, at least in principle, the route to hyperex- 
cit ability was discovere d on a molecular level. A recent migraine study by 



Vanmolkot et all (120071 ) showed a mutation in the gene coding the sodium 
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channel that leads to a slowed-down inactivation and a two-fold faster recov- 
ery from inactivation. In our simplified FHN system this is modeled by slowed 
inhibition, i.e., lower e values. As a consequence, the physiological state (A, 
Fig. 4d) of the cortex is shifted into hyperexcitable state (■). By definition, the 
physiological state ( A) is on a lower TAR- level than the hyperexcitable state 
(■). The therapeutic goal would be to shift the cortical state back to low a 
TAR-level, at which the pathological activity, i. e., the neurological symptoms, 
do not spread. 

To begin with, we need to specify dose response and toxicity curves. They 
are chosen arbitrarily to complete the example but without direct relation to 
the above mentioned subtype of migraine. Suppose that we have a drug with 
concentration acting exclusively on e (e. g., = c)- dose response graph 
is shown in Fig. 4b. It is similar to the one in Fig. 4a but rotated such that the 
response axis corresponds to the e-axis. Let another drug ('^ act exclusively 
on /3 (e. g., = /?). Its dose response graph is shown in Fig. 4c. It is mirrored 
(from Fig. 4a) such that the response axis corresponds to the /3-axis. For the 
sake of readability we replace the set {1,2} of indices i by {e, /?} and use 
for both the drug concentration and the drug name. Let the effective doses 50 
ECIq be the same for both (without loss of generality we choose ECI^q = 0.1), 
and let the toxic doses 50 TCIq differ. Usually TC^q is several times higher 
then the effective dose. We choose TQo = 5 and TCfo = 10. The drug C'^ 
has therefore a relatively higher therapeutic index TP = TCI^/ECIq than 
C^. Furthermore, we assume that both drugs have a sufficient efficacy, that is, 
they are each potent enough to drag the FHN system back to it physiological 
excitability defined by the TAR-isoline (Fig. 4d). The maximal response rates 
iimax) for both are chosen such that already 80% of is effective. Likewise 
and again without loss of generality we fix T^^^ at ^^ax- 

Given the example described above, we can calculate the cost for individual 
and combined therapy. If exclusively drug C,^ is administered to control ex- 
citability, the FHN system moves in parameter space along the /3-axis (down- 
wards in Fig. 4d) until the physiological TAR level (TAR low) is reached. The 
toxicity imposed by the therapeutic benefit in this case is 3.85% of the max- 
imal toxic level of . An exclusive administration of drug C*^ (which has a 
only half the therapeutic index TP compared to TI^) cost 7.4% of the max- 
imal toxic level T'^. This corresponds to a shift from point marked with ■ to 
the point marked with A in Fig. 4d. At last, we consider how the two drugs 
must be administered in a combined therapy. An optimal combined therapy 
reaches the physiological TAR level at 1.84% of the maximal combined toxic 
level Ttot = T'^ + T'^ (see the two black lines starting at (•) and terminating at 
the both effect axes in Fig. 4b and d). Due to the nonlinearity in the response 
and toxcity curves the total cost of the combined theraby is much less than 
the average of the individual therapies. 
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6 Discussion 



The phenotypical scheme of cortical spreading depression describes transient 
waves of massi ve depolariza t ion o f neurons and astrocytes. Such a spreading 
depolarization ( iDreier et all |2006| ) is associated with migraine, epilepsy and 
stroke. Their etiologies are mainly discussed in the context of hyperexcitability 
and disorders known as channelopathies, that is, diseases caused by a mutation 
in gene coding for ion channels. The challenge is to bridge the gap between 
the molecular level of the cause and macroscopic tissue level of the effects. 



We suggest that there are at least two independed routes towards a hyperex- 
citable state that supports transient wave forms in cortical tissue: one chang- 
ing the ratio of kinetic rates and one lowering the threshold. While the latter 
route changes the nuUclines, the former changes only the trajectories in the 
phase space. Consequently, there are also two routes out of the hyperexcitable 
regime. This is of particular interest when a critical therapeutic time window 
exits in which the volume of affected tissue is largely increased, as for exam- 
ple in peri-infarct depolarization. Then the therapeutic aim would be first to 
prevent tissue loss within the given therapeutic time window and re-establish 
the physiological value later. In Fig. 4d this would correspond to a path from 
the hyperexcitable state (■) directly to a state of low tissue at risk (•) and 
from there back to the physiological state (A). 



We have described efficient control of excitability by a simple pharmaco- 
dynamic model. In general, the study of beneficial effects of independent 
pathways is complicated by numerous interactions between pharmacokinet- 
ics, pharmacodynamics, and homeostatic factors and by individual variability. 
For example, the introduction of antagonistic behavior between and will 
complicate the geometrical structure of U. In general the coordinate system of 
the costs [T^, . . . ,T"'] will not build a Cartesian coordinate system. However, 
if we have a mapping from general costs and effects defined by Eqs. (5) and 
(7) we still can infer a metric structure of the parameter space of FHN. 



In conclusion, as a first step to implement the proposed control strategy ac- 
curate estimates given by pharmacokinetic-pharmacodynamic models of drug 
activity are needed for the regime of transient wave forms. We want to stress 
this point, namely that this estimates must be obtained in the regime of tran- 
sient wave forms, that is, a regime close to a bifurcation. Attempts to obtain 
estimates of pharmacokinetic-pharmacodynamic parameters in the regime far 
from this bifurcation will probably lead to inaccurate results. Then the next 
step is to construct an appropriate metric on the differential manifold of the 
parameter space in order to optimize the path from the hyperexcitable state 
back to the physiological state. 



10 



Acknowledgments 



We are grateful to G. Bordyugov, and V. Zykov for fruitful discussions. MAD 
would like to thank the Deutsche Forschungsgemeinschaft for financial support 
(DA-602/1-1). This work was also supported within the framework of Sfb 555. 



References 

Aurora, S. K., Ahmad, B. K., Welch, K. M., Bhardhwaj, P. & Ramadan, 
N. M. [1998] "Transcranial magnetic stimulation confirms hyperexcitability 
of occipital cortex in migraine," Neurology 50, 1111. 

Back, T. [1998] "Pathophysiology of the ischemic penumbra-revision of a con- 
cept," Cell Mol. NeuroUol. 18, 621. 

Baird, A. E., Benfield, A., Schlaug, G., Siewert, B., Lovblad, K. O., Edelman, 
R. R. & Warach, S. [1997] "Enlargement of human cerebral ischemic lesion 
volumes measured by diffusion-weighted magnetic resonance imaging," Ann. 
Neurol. 41, 581. 

Bures, J., Buresova, O. & Krivanck, J. [1974] The mechanism and applications 
of Ledo's Spreading Depression (Acadcmia, New York). 

Dahlem, M. A. & Chronicle, E. P. [2004] "A computational perspective on 
migraine aura," Prog. Neurobiol. 74(6), 351. 

Dahlem, M. A. & Miiller, S. C. [2004] "Reaction-diffusion waves in neuronal 
tissue and the window of cortical excitability," Ann. Phys. 13(7), 442. 

Dreier, J. P., Woitzik, J., Fabricius, M., Bhatia, R., Major, S., Drenckhahn, C., 
Lehmann, T.-N., Sarrafzadeh, A., Willumscn, L., Hartings, J. A., Sakowitz, 
O. W., Seemann, J. H. et al. [2006] "Delayed ischaemic neurological deficits 
after subarachnoid haemorrhage are associated with clusters of spreading 
depolarizations," Brain 129, 3224. 

Fabricius, M., Fuhr, S., Bhatia, R., Boutelle, M., Hashemi, P., Strong, A. J. 
& Lauritzcn, M. [2006] "Cortical spreading depression and peri-infarct de- 
polarization in acutely injured human cerebral cortex," Brain 129, 778. 

Fisher, M. & Takano, K. [1995] "The penumbra, therapeutic time window and 
acute ischaemic stroke," Baillieres Clin. Neurol. 4, 279. 

FitzHugh, R. [1961] "Impulses and physiological states in theoretical models 
of nerve membrane," Biophys. J. 1, 445. 

Giorgi, F. S., Lazzeri, G., Natale, G., ludice. A., Ruggieri, S., Paparelli, A., 
Murri, L. & Fornai, F. [2006] "MDMA and seizures: a dangerous liaison?," 
Ann. N. Y. Acad. Sci. 1074, 357. 

Gorji, A. & Spcckmann, E.- J. [2004] "Spreading depression enhances the spon- 
taneous epileptiform activity in human neocortical tissues," Eur. J. Neu- 
rosci. 19, 3371. 

Hadjikhani, N., Sanchez Del. Rio., M., Wu, O., Schwartz, D., Bakker, D., Fis- 



11 



chl, B., Kwong, K. K., Cutrer, F. M., Rosen, B. R., Tootell, R. B., Sorensen, 
A. G. & Moskowitz, M. A. [2001] "Mechanisms of migraine aura revealed 
by functional MRI in human visual cortex," Proc. Natl. Acad. Sci. USA 98, 
4687. 

Hakim, A. M. [1998] "Ischemic penumbra: the therapeutic window," Neurology 
51(3 Suppl 3), 44. 

Hakim, V. & Karma, A. [1999] "Theory of spiral wave dynamics in weakly 
excitable media: asymptotic reduction to a kinematic model and applica- 
tions," Phys. Rev. E 60, 5073. 

Hille, B. [1997] Ionic Channels of Excitable Membranes (Sinauer As). 

Hodgkin, A. L. & Huxley, A. F. [1952] "A Quantitative Description of Mem- 
brane Current and its Application to Conduction and Excitation in Nerve," 
J. Physiol. 117, 500. 

Hossmann, K. A. [1994] "Viability thresholds and the penumbra of focal is- 
chemia," Ann. Neurol. 36, 557. 

Kager, H., Wadman, W. J. & Somjen, G. G. [2000] "Simulated seizures and 
spreading depression in a neuron model incorporating interstitial space and 
ion concentrations," J. Neurophysiol. 84, 495. 

Kuznetsov, Y. A. [1995] Elements of Applied Bifurcation Theory (Springer, 
New York). 

Lashley, K. [1941] "Patterns of Cerebral Integration inicated by Scotomas of 
Migraine," Arch. Neurol. Psychiatry 46, 331. 

Lauritzen, M. [1994] "Pathophysiology of the migraine aura. The spreading 
depression theory [see comments]," Brain 117, 199. 

Leao, A. A. P. [1944] "Spreading depression of activity in the cerebral cortex," 
,/. Neurophysiol. 7, 359. 

Martins-Fcrrcira, H., Nedergaard, M. & Nicholson, C. [2000] "Perspectives on 
spreading depression," Brain. Res. Brain Res. Rev. 32(1), 215. 

Mayevsky, A., Doron, A., Manor, T., Meilin, S., Zarchin, N. & Ouaknine, G. E. 
[1996] "Cortical spreading depression recorded from the human brain using 
a multiparamctric monitoring system," Brain Res. 740, 268. 

Mihaliuk, E., Sakurai, T., Chirila, F. & Showalter, K. [2002] "Feedback stabi- 
lization of unstable propagating waves," Phys. Rev. E 65(6), 065602. 

Mikhailov, A. S. & Showalter, K. [2006] "Control of waves, patterns and tur- 
bulence in chemical systems," Phys. Rep. 425, 79. 

Mikhailov, A. S. & Zykov, V. S. [1991] "Kinematical theory of spiral waves 
in excitable media: comparison with numerical simulations," Physica D 52, 
379. 

Nagumo, J., Arimoto, S. & Yoshizawa., S. [1962] "An active pulse transmission 
hne simulating nerve axon.," Proc. IRE 50, 2061. 

Nedergaard, M. & Astrup, J. [1986] "Infarct rim: effect of hyperglycemia on di- 
rect current potential and [14C]2-deoxyglucose phosphorylation," J. Cereb. 
Blood Flow Metab. 6, 607. 

Reggia, J. A. & Montgomery, D. [1994] "Modeling cortical spreading depres- 
sion," Proc. Annu. Symp. Comput. Appl. Med. Care , 873. 



12 



Sakurai, T., Mihaliuk, E., Chirila, F. & Showalter, K. [2002] "Design and 
control of wave propagation patterns in excitable media," Science 296, 2009. 

Schindler, K. A., Leung, H., Lehnertz, K. & Elger, C. E. [2007] "How gener- 
alized are secondarily "generalized" tonic-clonic seizures?," J Neurol Neu- 
rosurg Psychiatry . 

Shapiro, B. E. [2001] "Osmotic forces and gap junctions in spreading depres- 
sion: a computational model," J. Comput. Neurosci. 10, 99. 

Somjen, G. G. [2001] "Mechanisms of spreading depression and hypoxic 
spreading depression-like depolarization," Physiol. Rev. 81, 1065. 

Somjen, G. G. & Aitken, P. G. [1984] "The ionic and metabohc responses 
associated with neuronal depression of Leao's type in cerebral cortex and 
in hippocampal formation," An. Acad. Bras. Cienc. 56(4), 495. 

Striimpel, C., Astrov, Y. A. & Purwins., H. G. [2002] "MuhiosciUatory pat- 
terns in a hybrid semiconductor gas-discharge system.," Phys. Rev. E 65, 
066210. 

Tobiasz, C. & Nicholson, C. [1982] "Tetrodotoxin resistant propagation and ex- 
tracellular sodium changes during spreading depression in rat cerebellum," 
Bram Res. 241, 329. 

Tuckwell, H. C. & Miura, R. M. [1978] "A mathematical model for spreading 
cortical depression," Biophysical J. 23, 257. 

Vanmolkot, K. R. J., Babini, E., de Vries, B., Stam, A. H., Frcilingcr, T., 
Terwindt, G. M., Norris, L., Haan, J., Frants, R. R., Ramadan, N. M., 
Ferrari, M. D., Pusch, M. et al. [2007] "The novel p.L1649Q mutation in 
the SCNIA epilepsy gene is associated with familial hemiplegic migraine: 
genetic and functional studies.," Hum. Mutat. 28, 522. 

Welch, K. M. [2005] "Brain hypcrcxcitability: the basis for antiepileptic drugs 
in migraine prevention," Headache 45, 25. 

Winfree, A. T. [1991] "Varieties of Spiral Wave Behaviour: An experimental- 
ist's approach to the theory of excitable media," Chaos 1, 303. 



13 




Fig. 1. The nullclines (bold) ii = and i; = in the phase space of the homogeneous 
FHN system with 7 = 0,/?= 1.4. Their intersection at {ug, Vg) is a stable fixed point. 
Three trajectories are drawn for e = 0.04: one canard trajectory (dotted), passing 
through the maximum of the nullisoclinc u = 0, and two trajectories starting at 
V = Vs nearby but on opposite sides of the canard trajectory. They diverge sharply, 
producing threshold behavior: (dashed) = super-threshold and (dashed-dotted) = 
sub-threshold stimulation. The parameter /? correlates with the threshold size, while 
A is in a certain range inversely related and therefore correlates with the excitability 
of the system (see text). 
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Fig. 2. Parameter space of the FHN system (7 = 0, (5 = 0). Besides the parameter 
e, a threshold parameter [5 and an excitability parameter A is used in (a) and (b), 
respectively, to span the parameter space (see Fig. 1). Three regimes exist defined 
by the spatio-temporal patterns that occur: waves, transient wave forms (TWF) 
and no spreading activity. The arrows mark two paths which are perpendicular to 
the border dS in (a) but not in (b). 
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Fig. 3. (a) Parameter space with tissue-at-risk(TAR)-isolines (dashed). The extent 
a TWF spreads is defined from the stimulation border to the location where the 
maximum activator concentration u hes above the nulhsochne -u = (inset). The 
red curve in the inset is the projection of a TWF from the infinite phase space 
of the space-dependent FHN system into the one of the homogeneous system for 
/3 = 1.4, e = 0.04. The projection is taken at the moment which we defined as 
the coUapse of the TWF. (b)-(d) Space-time plots of the transient wave forms 
(TWF) following a stimulation (increase of n by 2 for < x < 1, starting from 
the fixed point). Parameters: (b) e = 0.04, /3 = 1.4, (c) e = 0.05, /? = 1.4, (d) 
e = 0.095, /5 = 1.185. = 1 and 6 = for all parts. The color code denotes the 
activator u as defined in the inset of panel (a). 
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Fig. 4. (a) Dose response curves for an effective dose and a toxic dose witli tlie 
typical sigmoidal shape caused by tlie logaritlimic scaling. (b),(c) Same as in (a), 
but with rotated and inverted axes for two drugs acting on e and acting on /3. 
(d) Schematic parameter space with two TAR-isolines (TAR high and TAR low). A 
route to hyperexcitability can be caused by a shift in e (from ▲ to ■). The optimum 
path from the hyperexcitable state (■) towards the isoline of the physiological state 
(TAR low) using the two dose and toxic response curves in (b) and (c) leads to (•). 



17 



